Good views of the Galaxy 



Hiranya V. Peiris 

Princeton University Observatory, Peyton Hall, Princeton, NJ 08544; 
hiranya@astro.princeton.edu 

ABSTRACT 

Fitting Galactic structure models to star counts only provides useful 
information about the Galaxy in some directions. In this paper, we investigate 
the use of \ 2 goodness-of-fit tests to discriminate between degenerate Galactic 
structure models, and the implications of this technique for the Galactic spheroid 
and thick disk components. The axis ratio of the Galactic spheroid and the 
normalization of spheroid stars with respect to disk stars introduce a degenerate 
effect which means that Galactic structure models with certain combinations of 
these parameters are indistinguishable from each other in most directions. We 
present an analysis of the optimal directions in which these degeneracies can be 
lifted. Poisson and magnitude errors are taken into account, and an attempt 
is made to place an upper limit on the systematic error due to separation of 
spheroid stars from thick/old disk stars. We find that the magnitude range 
20 < V < 21 is the best for lifting most degeneracies, and present the optimal 
combinations of directions using which this can be achieved. We also give 
directions in which the signature of the presence of a Galactic thick disk can be 
most readily identified, and the directions in which contamination from a thick 
disk can be minimized. It is hoped that forthcoming data from large-scale sky 
surveys would reveal much about the structure of our Galaxy using star count 
techniques. 



Subject headings: Galaxy: structure — stars: statistics 



1. Introduction 

Since we observe our Galaxy from within it, we must use indirect tools such as star 
counts to probe its structure. The success of this method relies on assuming that the 
components of the Milky Way Galaxy have stars distributed similarly to those of galaxies of 
the same Hubble type, and comparing models with observations only in regions where little 
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obscuration is known to be present. Observed "star counts" (i.e. the distribution of number 
counts of stars in apparent magnitude and color within a given area/direction in the sky) 
are then used to determine parameters in distribution functions whose overall shapes are 
assumed known. 

Within the past two decades, Galactic structure models of varying degrees of complexity 
have been developed (Bahcall & Soneira, 1980, 1981, 1984, Gilmore 1984] , Robin and Creze 



1986a, 1986b, |Reid fc Majewski 1993| ). However, during the same period, relatively few sets 



of observational star count data have been published. 

Galactic structure models, for a given set of structure parameters, typically predict 
differential number counts (stars / apparent magnitude bin / area) and frequency- color 
distributions (stars / color bin / area) for a specified direction. Because the Galaxy consists 
of at least two components, obeying different density laws and being characterized by 
different stellar populations, these distributions change with the color range, magnitude 
range and direction under consideration. It is important to match both the shapes of these 
distributions and the absolute numbers of stars with observations. 

Existing models predict differential number counts reasonably well, to faint magnitudes 
of V ~ 21, but recent studies suggest that agreement between observed and predicted 
frequency-color distributions breaks down at faint magnitudes below V ~ 19 ( |Rcid| 
|fc Majewski 19931) . "Standard" Galactic structure models appear to overestimate the 
contribution of (blue) spheroid stars substantially, and the contribution of (red) disk stars 
is correspondingly underestimated in order to match the total number counts. Contrary to 
prediction, observations show that disk stars are the majority population to at least V ~ 21. 
However, the small sizes of fields in which data are available mean that comparisons in 
frequency-color space are complicated by uncertainties due to small-number statistics. 

A further hindrance to the determination of distribution function parameters arises 
due to the fact that data are only available in a small number of directions. This presents 
a problem in trying to fit model parameters which cause degenerate effects in the number 
counts. For instance, increasing the normalization of the Galactic spheroid and decreasing 
its axis ratio b / a (or vice versa) represents a degeneracy. This means that in most 
directions, one cannot distinguish between models with a flattened spheroid plus a low 
normalization ratio of spheroid stars to disk stars, and those with a high axis ratio plus a 
high normalization. 

We anticipate the availability of star count data from large-scale, deep sky surveys in 
the near future (see § [| and § ||) which would solve the problem of excessive statistical 
noise, and provide data in almost any direction we desire. In this paper, we investigate the 
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possibility of tackling the degeneracy problem by attempting to find two or more directions 
in which x 2 goodness-of-fit tests can distinguish between spheroid parameter combinations 
that yield degenerate results in general. 

A degenerate situation is defined as one in which, in a general direction, there is 
no significant difference in results given by a model where a pair of spheroid parameters 
are characterized by values (p, q) and a model characterized by values (r, s) . This means 
that one cannot learn very much about these spheroid parameters even if these models 
fit the spheroidal star counts in that direction with reasonable accuracy. The degenerate 
parameters we examine here are the spheroid axis ratio and the normalization of spheroid 
stars to disk stars in the solar neighborhood. 



By using x 2 tests with a specified set of directions as "bins" (see e.g. Press et al. 1992| ) 



we have investigated the optimal directions for breaking this degeneracy. These tests are 
based on fitting the total number of spheroid stars in each direction "bin" , for a given set of 
directions and a particular magnitude range. Therefore, both the observed number count 
and color information is used, since the spheroid counts have to be separated from the total 
counts using the frequency-color distribution, as explained in § |4 . 1 . 4j . 



§ H describes the Galactic structure models we have used in these tests. Details of the 
directions used in the analysis are found in § j|. A description of the x 2 tests, along with 
our results, are presented in § |j. We conclude, in § |5], that there are two or more directions 
in which most "degenerate" model pairs can be distinguished. We find the magnitude range 
20 < V < 21 to be the most useful for this purpose. We also make suggestions for the 
best directions in which to separate spheroidal stars from old- and thick-disk stars using a 
frequency-color distribution diagram, and conversely, the directions in which the presence 
of a thick disk would be most clearly observed. 



2. The Models 

While it is our intention to illustrate the use of this technique in general, for specificity 
we will use the model that is usually referred to in the literature as the Bahcall-Soneira 
(B&S) model. We have taken the realization of the B&S model in the form of the Export 
Code that is available on the Internet ( [http: / /www.sns.ias.edu/^jnb/Html/galaxy.html| ). 
This model is described in Bahcall (1986). We have added a few refinements as described 
below. 

The disk luminosity function (LF) has been modified between 9.5 < My < 18.0 to 
match the LF recently derived from Hubble Space Telescope (HST) M dwarfs (Figure 2, 
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Gould et al. 199? ). A "sech 2 " model has been used for the vertical density distribution 



function of main sequence disk stars in the same absolute magnitude range, as given in Eq. 
3.2 of the same work. The HST LF was derived using two different methods: a maximum 
likelihood (ML) fit that takes into account both measurement errors and Malmquist bias, 
and a naive binning method where the total number of stars in a given magnitude bin is 
divided by the effective volume integrated over that bin. This binning method takes into 
account neither Malmquist bias nor observational error. The two alternative derivations 
agree quite well over most of the LF, but at the faint end (13.5 < My < 18.5) the 
ML procedure could by affected by a statistical fluctuation. Because Poisson errors are 
potentially a much more serious problem at the faint end than Malmquist bias (since the 
HST survey extends to the "top" of the disk), we choose to use the form of the LF derived 
using the binning method. 

The spheroid LF has been modified between 7.5 < My < 13.5 to match another LF 
derived from HST star counts (Figure 3, pould et al. 1998]) . Note that this LF differs in its 



zero point by a factor of ~ 2 with studies of kinematically-selected local spheroid subdwarfs 
( Palm et al. 1995|) . 



The disk and spheroid LFs used are shown in Figure [I]. In addition to the spheroidal 
parameters whose values are described in § |4], the color-magnitude diagram (CMD) of M13 
was adopted for the spheroid CMD. 

In addition to this two-component model, a three-component model was created with 
the addition of a thick disk with a scale height of 1200 pc, a 47 Tuc-like CMD, a disk-like 
LF, and a normalization to disk counts at the solar neighborhood of 2%. This was used in 
analyzing the separation error in identifying spheroid counts from disk/thick-disk counts 
using frequency-color diagrams, as described in § |4.1.4| . 

The density laws used for the old disk, the spheroid and the thick disk are summarized 
in Table |], which also gives the normalizations used. 



3. The Directions 

We selected a representative set of directions in which to carry out the analysis, which 
are shown in Table 0. 

Directions 1 & 2 were selected following a quantitative analysis of directions in 
the (I 11 = 0° -> 180°) and (l n = 90° -> 270°) planes, as the optimal directions in 
which to distinguish between the effects of the normalization and the axis ratio of the 
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spheroid. Bahcall and Soneira (1980) explain why these planes are especially important 
for determining the characteristics of the spheroid. We ran the model at 5° intervals in 
b 11 in each of these planes, and calculated the variation with b 11 of the ratio of spheroid 
counts to disk counts, and the ratio of spheroid counts for pairs of near-degenerate models: 
for example, using the notation (disk : spheroid normalization ratio, spheroid axis ratio), 
(800 : 1, 1.0) : (500 : 1,0.8). We then looked for directions with a high spheroid : disk ratio 
(in order to minimize the Poisson error in total spheroid counts as well as the error due 
to contamination of spheroid counts by disk stars) while at the same time, maximizing 
the quantity | (800 : 1,1.0) : (500 : 1,0.8) - 1 |. Directions 1 & 2 optimized both these 
quantities. 

Directions 3-10 were selected because fields in these directions are being prepared by 
the Digital Palomar Observatory Sky Survey (DPOSSII; S. G. Djorgovski & S. C. Odewahn 
1998, private communication; Djorgovski et al. 1997, 1998) for star count studies. It would 
then be possible to see how the techniques discussed here in a theoretical sense work in 
practice. 



4. x 2 Analysis 



The x 2 statistic is 



2 J^, [A model (m u m 2 ,l,b)dQ, - A observed (m 1 ,m 2 ,l,b)dtt} 2 
X = 2^ -2 , I 1 ) 



dir=l 



where dir = (/, b) is a given direction, A(mi, 1712, 1, b)dQ is the total number of spheroid 
stars with apparent magnitudes in the range m\ < m < m 2 in the direction (/, b) per 
projected sky area dQ, and a is the error, which can have several components as described 
below. Reduced x 2 is given by x 2 / u where v is the number of degrees of freedom, in 
this case equal to the number of directions N (i.e. the total predicted counts are not 
renormalized to the same total counts as the data). A large value of reduced x 2 indicates 
that the null hypothesis that the A observe dS are drawn from the same population as the 
A mo deiS is rather unlikely. 

Since we do not currently have observed number counts in all the directions, we picked 
one model out of the set covering the investigated parameter-space to be the "observed" 
number counts and compared it against each of the rest of the models in the set (which 
represented the "theory"). The process was then repeated till each of the models had been 
picked as the one representing actual data. 
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We define the criterion for lifting the degeneracy of spheroid parameters as obtaining 
the probability of finding a value of reduced x 2 greater than or equal to the "observed" value 
< 10~ 5 . That is, we aim to reject the null hypothesis (that the "observed" and "theoretical" 
models are the same) at a confidence level of 0.001%. The range of models we used spanned 
the parameter space (disk : spheroid normalization ratio) x (spheroid axis ratio) = (500 : 
1,800 : 1) x (1.0,0.8,0.6,0.4). Thus there are 56 possible pairings of models picked as 
"observed" and "theoretical" . 

X 2 tests were carried out using the set of directions detailed in § |3| as "bins." The 
directions yielding the highest contribution to the x 2 value were noted as the optimal 
directions to use for breaking the degeneracy. It was then checked whether x 2 tests using 
just those optimal directions as bins would cause the degeneracy to be lifted. 

The error term a can have contributions from Poisson error, magnitude error, and 
the systematic "separation error" which arises because spheroid stars cannot be separated 
cleanly from old/thick disk stars in observations, thereby introducing an error into the 
total spheroid counts. In general, a magnitude error of order ~ 0.1 ma9 contributes less to 
the error budget than the Poisson error. In a "good" direction, the separation error should 
contribute the least to the error budget and the number counts are Poisson-noise-limited. 
In a "bad" direction, the separation error can be by far the greatest source of uncertainty. 



In § fLl| , we first consider the case where the spheroid counts can be cleanly separated, and 



then in § 4.2 we consider the effect of adding the separation error. 



4.1. Analysis excluding separation error 

X 2 tests were first carried out using the two-component (spheroid + old disk) model, 
taking into account the Poisson error of the counts and a magnitude error of 0.1 ma9 , 
but assuming that spheroid stars can be separated cleanly from old/thick disk stars in 
observations. 

The magnitude error was calculated as follows. Let N(m\ < m < m 2 ) represent the 
total spheroidal number counts for apparent magnitude range m\ < m < 7712, and let the 
error in the zero-point be 5m. Then we calculate N + (m\ + 5m < m < m 2 + 5m) and 
AL (mi — 5m < m < m 2 — 5m) . The error in the total spheroidal number counts due to the 
zero-point error in magnitude is then given by Max(\ N — N + \, \ N — AL |). 
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4-1.1. Directions 1 & 2 

Directions 1 and 2 in Table Q were found to be the best overall for discriminating 
between ordinarily degenerate models using y 2 tests. 

Using the notation [{model taken as observation), {model taken as prediction)], 
almost all model pairings were distinguishable from each other in the apparent 
magnitude range 19 < V < 20. The exceptions were [(500 : 1,0.8), (800 : 1,1.0)] and 
[(800 : 1,1.0), (500 : 1,0.8)]. Upon examining the range 20 < V < 21 the degeneracy 
of these model pairings were also lifted, while the other model pairings also remained 
non-degenerate. 



4.1.2. Directions 3-10 (19 <V < 20) 

Considering the directions 3-10 which are being prepared by the DPOSSII survey, 
X 2 tests in directions 3 and 4 in apparent magnitude range 19 < V < 20 can distinguish 
between all model pairings with the exceptions shown in Table ^. 



4.1.3. Directions 3-10 {20 <V <21) 

Considering the model pairings in the apparent magnitude range 20 < V < 21 in the 
DPOSSII directions, \ 2 tests can distinguish between all model pairings in directions 3 and 
4 except for those shown in Table |j. 



4-1-4- Discussion 

Observationally, spheroid stars can be separated from old/thick disk stars using the 
bi-modal distribution that appears in star color-frequency profiles in certain directions 
at faint magnitude ranges. The blue and red peaks consist of spheroid and disk stars, 
respectively. This distribution has been interpreted as arising because the disk and spheroid 
components have different density gradients in a magnitude-limited survey (|Bahcall| 



|fc Soneira 1980Q . The sharp density gradient in the disk for directions far from the 
Galactic plane favors relatively nearby and intrinsically faint (red) stars at faint apparent 
magnitudes. Conversely, the shallow density gradient in the spheroid favors relatively 
distant, intrinsically bright (blue) stars at the same faint apparent magnitudes, because the 
effective volume increases at large distances. 
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18 < V < 22 is the apparent magnitude range where this double-peaked distribution 
is most pronounced, enabling the cleanest separation of spheroidal stars in observations 
( |Bahcall fc Soneira 1980| , [Reid fc Majewski 1993| ). In this interval, the spheroid counts peak 



at B — V ~ 0.5, with a narrow range in color. This corresponds roughly to the main-sequence 
turn-off, at My ~ 4.5. For a given apparent magnitude interval m\ < m < rri2, only stars 
within a range of absolute magnitudes Mi < M < M 2 contribute to the observed number 
counts. The majority contribution to the blue spheroidal peak comes from stars near the 
main sequence turn-off, with +4 < My + 6 (corresponding to ~ 0.5 — > 1.0 M Q ). Evolved 
stars with My < 3.5 make negligible contribution to deep star counts at this magnitude 
range. 

The shallow density gradient in the spheroid leads to a broad distribution of number 
counts with distance in a given direction, peaking at ~ 7 kpc in the direction of the North 
Galactic Pole ( [Bahcall &: Soncira 1980 ). Therefore, the fainter the magnitude range under 



consideration, the better the performance of the \ 2 statistic, because spheroid number 
counts are still rising at V ~ 22. Thus the Poisson errors and magnitude errors at fainter 
magnitudes are smaller. This is the reason that the magnitude range 20 < V < 21 is more 
useful for lifting degeneracies using the \ 2 tests than the brighter range 19 < V < 20. Once 
we go even fainter to the range 21 < V < 22, photometric uncertainties and observational 
errors, such as star-galaxy separation and contamination of the blue spheroidal population 
by quasars and compact emission line galaxies (CELGs) increasingly start to affect counts. 
Thus the intermediate apparent magnitude range 20 < V < 21 is the optimal one to use for 
the purpose of lifting degeneracies by this method. 



4.2. Analysis including separation error 

Figures |2] and ^ show frequency distributions of star colors for a particular "standard" 
model, in the ten directions given in Table |2|, for magnitude ranges 19 < V < 20 and 
20 < V < 21. The maximum error in separating out the spheroid stars arises if there is a 
thick disk component which tends to fill up the valley between the double-peaks due to 
spheroid and old disk stars, or in the worst cases, peaks within the blue spheroid peak. The 
separation error is calculated as the area under the blue peak in total (spheroid + old disk 
+ thick disk) counts, minus the calculated (i.e. model) spheroid counts. The "blue peak" is 
taken to be bounded on the red side by the color corresponding to the minimum point of 
the total counts. 

Now, we consider the effect of including the systematic error due to separating spheroid 
stars from old/thick disk stars in observations, using the three-component (spheroid + old 
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disk + thick disk) model. The separation error is an overestimate since this contamination 
in observations can be corrected for to some extent by subtracting model disk counts from 
the total counts. 

The contamination due to the thick disk was found to dominate that due to old disk 
stars in all cases. This error was found to be smaller in the apparent magnitude range 
20 < V < 21 than in the range 19 < V < 20. Hence the following analysis was only carried 
out for the fainter range. Further, since the ratio of separation error to total spheroid 
counts increases significantly for models with smaller axis ratios, the models with axis 
ratio 0.4 were excluded from the analysis. Thus, model parameter space was restricted to 
(500 : 1,800 : 1) x (1.0,0.8,0.6). 

Upon including the separation error, directions 1 and 2, found to be the best overall 
in discriminating between model parameters, were found to be virtually useless, since the 
thick disk counts peaked within the spheroid peak for these directions. 

In the case of the DPOSSII directions (3-10), all the model pairings which could be 
distinguished in this apparent magnitude range using \ 2 tests in directions 3 and 4 remained 
non-degenerate, except [(800 : 1,0.6), (500 : 1,0.6)], which became degenerate. The rest 
of the model pairings which had been non-degenerate in directions given in Table |] also 
became degenerate. 

Table |^ shows whether these degenerate models can be discriminated using other 
combinations of directions. 

During this analysis, it was found directions 7 and 9 are generally good for minimizing 
the effects due to the thick disk. Conversely, directions 1 and 10 are particularly good for 
looking for the presence of the thick disk. 

5. Conclusions 

We found that directions 1 (I 11 = 0°, b 11 = 40°) and 2 (I 11 = 90°, b 11 = 40°) in apparent 
magnitude range 20 < V < 21 were the most effective in distinguishing between effects due 
to degenerate parameters, assuming that spheroid stars could be separated cleanly from 
disk and thick-disk stars. 

In the case that the maximum possible separation error (i.e. uncorrected for disk and 
thick disk counts) was included in the analysis, directions 3 (the North Galactic Pole) 
and 4 (I 11 = 67°, b 11 = 49°) in apparent magnitude range 20 < V < 21 were the most 
useful. In those models where these directions failed to lift the degeneracy, directions 
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7 (I 11 = 111°, b 11 = -46°) and 9 (I 11 = 172°, b 11 = 48°) were the best for minimizing 
contamination from the thick disk, which dominates the separation error. 

Directions 1 (I 11 = 0°,b n = 40°) and 10 (l n = Ql°,b n = -37°) were found to be the 
most effective in detecting the presence of a thick disk. 

This is an exciting time in which deep, high-quality data-sets covering large areas of 
the sky are about to become available from much-anticipated surveys such as the Sloan 
Digital Sky Survey (SDSS, see e.g. |Gunn fc Knapp 1993| , punn fc Weinberg 1995| ) and 
the DPOSSII (see § ^ for references). The use of survey data for star count studies is 
complicated by the fact that many existing "standard" Galactic structure models are 
based on color-magnitude diagrams (CMDs) and luminosity functions (LFs) determined 
in the standard Johnson-Morgan- Cousins UBV photometric system, and are only useful 
for comparison with data in the visual band or in photometric systems to which accurate 
transformations exist. However, in cases such as that of the SDSS, which uses the 
non-standard u'g'r'i' z' photometric system (Fukugita et al. 1996), such models will soon be 



re-created using LFs and CMDs determined in the SDSS photometric bands. The DPOSSII 
uses three photographic JFN bands calibrated to the Gunn gri bands. Though these 
are different from the photoelectric/CCD gri bands (Palomar, 4-Shooter, SDSS) they are 
well-defined. Once the work of creating a model based on the SDSS bands is completed, it 
should be fairly straightforward to translate the stellar sequences and LFs to the DPOSSII 
gri bands. Until such models are available, a different approach to star count models, such 
as the evolutionary stellar population synthesis technique ( Pan 1999] ) can be used. 

The possibility of fitting models simultaneously in a multitude of directions with large 
samples containing number, magnitude and color information should revolutionize studies 
of Galactic structure using star counts. 



The author is very grateful to John Bahcall for his helpful advice. She wishes to thank 
Andy Gould for providing machine-readable versions of Figure 2 of Gould et al. (1997) and 
Figure 3 of Gould et al. (1998); David Spergel and Xiaohui Fan for instructive discussions 
on Galactic Structure models; and Mark Jackson for useful comments on the manuscript. 



Table 1. The Assumed Stellar Distributions. 



Component 



Distribution 



Disk no = ud{Ro) exp[— z/H (My)] exp[— (x — Ro)/h] 

Spheroid n S ph = n sp h 

(i?o)(i?/i?o) _7/8 exp[-10.093( J R/i? ) 1/4 + 10.093] 
x 125(^^0) 6/8 cxp[-10.093(i?/i? ) _1/4 + 10.093], R < 0.03i? 
x[l - 0.08669/(i?/i? ) 1/4 ], R>0mR o 



Thick Disk utd = ^td(-Ro) ex p[— z/Htd] exp[— (x — Ro)/h] 

Normalization n D (R ) = 0.13pc" 3 , n TD (R ) = 0.0026pc~ 3 , 

n sp h{Ro) = 0.00026pc" 3 or 0.0001625pc~ 3 depending on model 



Note. — Here z is the distance perpendicular to the plane, x is the galactocentric 
distance in the plane, and h is the old disk scale length. Galactocentric distance 
R = (x 2 + z 2 /k 2 ) 1 / 2 , where k is the axis ratio and 1 — n is the ellipticity. We adopt 
Ro = 8 kpc and h = 3.5 kpc. The old disk scale height H(My) is given in Bahcall 
and Soneira (1980). The thick disk is taken to have a scale height Htd = 1-2 kpc 
and a 47-Tuc-like CMD. 
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Table 2. A Representative Set of Directions. 



Direction 


Z 7/ (1950) 


6 /7 (1950) 


Comments a 


1 


0° 


+40° 


BS9 b 


2 


90 


+40 


Between BS4 and BS7 b 


3 




+90 


Galactic pole (SA57, BS1) C <' 


4 


67 


+49 


BS16 c < d 


5 


180 


+50 


BS6 C 


6 


180 


+30 


BS8 C 


7 


111 


-46 


BS11, SA68 C 


8 


167 


-51 


BS13 C 


9 


172 


+48 


BS14 C 


10 


61 


-37 


BS17 C 



a The field identification BSn denotes the designation assigned to a given 
direction in Table 1 of Bahcall and Soneira (1981), where quantities that 
can be efficiently studied in that direction are listed. The Selected Area 
number is also given for fields which have been studied. 

b These directions were determined to be the best for distinguishing 
between effects due to spheroid flattening and normalization using the two- 
component (old disk + spheroid) standard model. 

c These fields are being made available by the DPOSSII survey for star 
count studies. 

d These directions were determined to be the best for determining the 
spheroid parameters out of the DPOSSII directions. 
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Table 3. Model pairings degenerate in directions 3 and 4 for 19 < V < 20 excluding 

separation error. 



Model pairing a Min. directions required to lift degeneracy 



[(500 : 


: 1,0.8), (800 : 


1,1.0)] 


6, 7 


[(800 : 


: 1,1.0), (500 : 


1,0.8)] 


6, 7, 10 


[(500 : 


: 1,0.6), (800 : 


1,0.8)] 


b 


[(800 : 


: 1,0.8), (500 : 


1,0.6)] 


b 


[(500 : 


: 1,0.4), (800 : 


1,0.6)] 


b 


[(800 : 


: 1,0.4), (500 : 


1,0.4)] 


4, 10 



a Uses notation [(model taken as observation), (model taken as prediction)]. 

b This model pairing cannot be distinguished using the DPOSSII directions 
3-10 in the apparent magnitude range 19 < V < 20. 



-14- 



Table 4. Model pairings degenerate in directions 3 and 4 for 20 < V < 21 excluding 

separation error. 



Model pairing 3, Min. directions required to lift degeneracy 

[(500 : 1,0.8), (800 : 1,1.0)] 6, 10 

[(800 : 1,1.0), (500 : 1,0.8)] 6, 10 

[(500 : 1,0.6), (800 : 1,0.8)] 3, 6 

[(800 : 1,0.8), (500 : 1,0.6)] ••• b 



a Uses notation [(model taken as observation), (model taken as prediction)]. 

b This model pairing cannot be distinguished using the DPOSSII directions 
3-10 in the apparent magnitude range 20 < V < 21. 
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Table 5. Model pairings degenerate in directions 3 and 4 for 20 < V < 21 including 

separation error. 



Model pairing a Min. directions required to lift degeneracy 



[(500 : 


1,0.8), (800: 


1,1.0)] 


7, 9 


[(800 : 


1,1.0), (500: 


1,0.8)] 


... b 


[(500 : 


1,0.6), (800: 


1,0.8)] 


b 


[(800 : 


1,0.8), (500: 


1,0.6)] 


b 


[(800 : 


1,0.6), (500: 


1,0.6)] 


7, 9 



a Uses notation [(model taken as observation), (model taken as prediction)]. 

b This model pairing cannot be distinguished using the DPOSSII directions 
3-10 in the apparent magnitude range 20 < V < 21 using x 2 tests including 
the separation error. 
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Fig. I. — The adopted disk and spheroid luminosity functions (LFs). The LFs adopted are 
the same form those described in Bahcall (1986) except for the following refinements: the 
disk LF for My > 9.5 has been modified to match smoothly the LF derived from HST star 
counts given in Figure 2 of Gould et al. (1997) and the spheroid LF for My > 7.5 has been 
modified to match the LF given in Figure 3 of Gould et al. (1998), also derived from HST 
data. 
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Fig. 2. — Frequency distribution of star colors for a "standard" model (500 : 1, 0.8) flBahcall 
1986|) plus a thick-disk as described in § ||, in the apparent magnitude range 19 < V < 20. 
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Fig. 3. — Frequency distribution of star colors for a "standard" model (500 : 1, 0.8) flBahcall 
1986|) plus a thick-disk as described in § ||, in the apparent magnitude range 20 < V < 21. 



